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Statement  of  Work 


Under  Contract  No.  19628-80-C-0020  D.P.  Kennedy  &  Assoc.,  Inc. 
has  undertaken  investigations  into  the  current-voltage  relationship  of  low 
barrier  Schottky  diodes.  A  particular  objective  of  this  program  was  to  de¬ 
velop  understanding  of  the  large  reverse  current  observed  in  Pt-Si  diodes 
on  p-type  silicon,  when  operating  at  low  temperatures  (77°  to  110°  k). 
Additionally,  after  identification  of  appropriate  physical  mechanisms  it 
was  desired  to  undertaken  rigorous  numerical  investigations  of  the  topic. 
Such  numerical  investigations  could  offer  a  means  for  quantitative  com¬ 
parisons  between  theory  and  measurements  upon  experimental  Pt-Si  semi¬ 
conductor  devices . 
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Summary  and  Malor  Results 


The  topic  of  this  investigation  was  based  upon  experimental  ob¬ 
servations  of  an  exceedingly  large  reverse  electric  current,  as  exhibited 
by  p-type  Pt-Si  Schottky  barrier  diodes  when  operating  at  a  low  ambient 
temperature  (77  °-  110°K).  Numerous  explanations  of  this  phenomenon 
are  based  upon  a  proposed  mechanism  of  hole  generation  within  metal¬ 
lurgical  regions  of  this  device.  It  has  been  suggested  that  such  holes 
transverse  the  interface  barrier,  much  like  conduction  electrons  in  an 
n-type  Schottky  barrier  diode  and,  thereafter,  produce  the  observed 
electric  current.  Because  hole  generation,  in  the  classical  sense,  can¬ 
not  exist  in  a  metallic  conductor,  this  investigation  was  directed  toward 
other  mechanisms;  mechanisms  capable  of  producing  a  large  reverse  elec¬ 
tric  current. 

Calculations  show  that  an  insignificant  number  of  mobile  electrons 
exist  in  p-type  silicon  at  low  ambient  temperatures  (77°  to  110° K)  and,  there 
fore,  the  observed  reverse  electric  current  is  unlikely  to  arise  from  normal 
hole/electron  pair  generation  within  the  silicon  space-charge  layer.  As 
a  consequence,  the  observed  electric  current  is  suggested  to  arise  from 
mobile  hole  generation  at  the  silicon-metal  interface  through  mechanisms 
heretofore  given  little  recognition  in  the  technical  literature.  Specifically, 
it  is  proposed  that  silicon  valence  electrons  enter  unoccupied  electron 
states  within  the  metallic  conductor  and,  thereby,  produce  mobile  holes 
without  a  traditional  hole-electron  pair  generation  in  the  semiconductor 
material. 

A  study  of  the  technical  literature  was  undertaken  on  the  topic 
of  electron  field-emission  from  silicon.  Some  experimental  data  on  elec¬ 
tron  field  emission  has  been  interpreted  as  arising  from  electron  emis¬ 
sion  from  the  valence  band  of  p-type  silicon.  Additionally,  using  ap¬ 
proximation  methods  other  workers  have  shown  that  field  emission  from 
the  valence  band  should,  indeed,  be  expected  from  p-type  silicon. 
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It  is  important  to  note  that  these  studies  of  field-emission  from  silicon 
were  undertaken  with  no  apparent  recognition  of  applicable  mechanisms 
for  p-type  Schottky  barrier  diodes . 

From  traditional  theory  of  metallic  conductors,  it  is  shown  that 
a  large  number  of  unoccupied  energy  states  exist  below  the  Fermi-level. 
Additionally,  it  is  shown  that  interface  valence  electrons  can  enter  these 
unoccupied  states  with  a  degree  of  excitation  below  that  required  by  clas¬ 
sical  theory  for  this  type  of  electron  transition.  As  a  consequence,  an 
effective  energy  barrier  could  be  defined  for  a  p-type  Schottky  barrier 
that  is  less  than  that  given  by  traditional  theory  of  metal- semiconductor 
contacts . 

For  a  p-type  Schottky  barrier,  detailed  calculations  are  presented 
for  the  reverse  electric  current.  These  calculations  are  based  upon  tradi¬ 
tional  concepts  drawn  from  the  statistical  theory  of  matter.  Specifically, 
the  rate  is  calculated  at  which  electron  transitions  take  place  from  the 
semiconductor  valence  band  to  the  conduction  band  of  the  metallic  con¬ 
tact.  Such  calculations  have  been  undertaken  for  a  large  assumed  range 
of  conduction  band  free  electron  density;  a  range  that  would  include  both 
degenerate  semiconductor  material  and  an  ideal  metal. 

Resulting  from  these  calculations  is  much  understanding  about  the 
mechanisms  involved  in  the  operation  of  a  p-type  Schottky  barrier.  It  is 
demonstrated  that  despite  the  experimental  observation  of  bidirectional 
electric  current,  a  true  interface  energy  barrier  exists  in  only  a  limited 
type  of  p-type  semiconductor  device.  In  particular  a  small  reverse  elec¬ 
tric  current  can  arise  from  one,  or  both,  of  two  different  mechanisms. 

A  small  reverse  current  can  arise  from  the  distribution  of  unoccupied  energy 
states  in  the  metal  conduction  band,  relative  to  the  energy  band  structure 
of  the  p-type  semiconductor  material.  Additionally,  an  energy  barrier  in 
the  classical  sense  can  be  shown  when  the  conduction  band  edge  of  the 
metal  is  at  a  higher  energy  than  the  valence  band  edge  of  the  semicon¬ 
ductor.  This  second  mechanism  will  probably  exist  only  if  the  adjacent 


material  is  a  highly  degenerate  semiconductor,  rather  than  a  metal.  In 
line  with  the  basic  objectives  of  this  research  ,  many  calculations  were 
undertaken  on  the  expected  reverse  current  from  low  barrier  semiconduc¬ 
tor  devices - defining  the  interface  energy  barrier  in  a  manner  consis- 

tant  with  current  textbook  concepts.  It  is  shown  that  relatively  poor 
rectification  characteristics  should  exist  for  low  barrier  devices,  if 
the  contacting  material  is  a  true  metal.  If,  instead,  the  contacting 
material  is  a  degenerate  n-type  semiconductor  the  rectification  prop¬ 
erties  could  be  substantially  improved,  although  in  most  cases  the  expected 
reverse  current  will  be  larger  than  calculated  from  classical  concepts  of 
Schottky  barrier  operation . 

For  a  rigorous  study  of  this  topic  a  computer  program  was  de¬ 
veloped  that  is  capable  of  calculating  the  electrostatic  potential  in  a 
reverse  biased  Schottky  barrier  diode,  in  the  presence  of  a  large  elec¬ 
tric  current.  This  computer  program  involves  the  simultaneous  solution 
of  Poisson's  equation  and  the  mobile  hole  continuity  equation.  Two 
distinctly  different  hole  generation  mechanisms  are  utilized:  first,  the 
classical  HSR  hole-electron  generation  through  a  single  mid-gap  recomb¬ 
ination  center  and,  second,  hole  generation  from  the  top  of  the  valence 
band.  This  computer  program  is  capable  of  accurately  calculating  the 
Schottky  barrier  potential  distribution,  even  in  the  presence  of  a  large 
reverse  electric  current .  Insufficient  time  was  available  during  the 
present  contract  period  to  complete  coupling  this  rigorous  solution 
to  the  calculation  of  electron  transition  at  the  metal/semiconductor  in¬ 
terface  . 
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Chapter  I 


Theory  of  Operation  for  a  p-type  Schottky  Barrier 

1.0  Introduction 

Experiment  shows  that  low-barrier  Schottky  barrier  diodes  on  p-type 
silicon  exhibit  poor  rectification  qualities,  even  at  low  ambient  tempera¬ 
tures.  Such  experiments  have  directed  workers  toward  a  rigorous  study 
of  physical  mechanisms  associated  with  this  type  of  semiconductor  device. 

In  particular,  the  physical  mechanisms  contributing  to  a  large  reverse  elec¬ 
tric  current.  Therein  we  have  the  basis  and  direction  of  the  present  re¬ 
search  on  Schottky  barrier  operation.  A  principle  goal  for  this  research  is 
to  undertake  a  mathematical  study  of  the  physics  associated  with  the  poor 
reverse  current  characteristics  of  the  Schottky  barrier  diode,  placing  particu¬ 
lar  emphasis  upon  devices  composed  of  Pt-Si  on  p-type  silicon. 

An  initial  survey  of  the  technical  literature  shows  that  most  theory 
of  Schottky  barrier  diode  operation  is  based  upon  devices  containing  n-type 
semiconductor  material.  This  theory  is  in  detail  and  mathematically  rig¬ 
orous.  Furthermore,  comparisons  between  theory  and  experiment  are  in 
qualitative  agreement,  and  many  are  in  quantitative  agreement,  for  most 
of  the  important  characteristics  exhibited  by  this  device.  A  source  of 
difficulty  in  such  comparisons  has  always  been  a  lack  of  detailed  know¬ 
ledge  about  the  structure  of  the  metal-semiconductor  interface.  Theoreti¬ 
cal  investigations  are  frequently  based  upon  an  idealization  of  this  inter¬ 
face,  an  idealization  that  renders  the  problem  tractable  from  a  mathemati¬ 
cal  and  physical  point  of  view. 
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Nonetheless,  it  Is  believed  safe  to  say  that  an  extensive  under¬ 
standing  is  available  on  the  operation  of  n-type  metal-semiconductor 
contacts  of  the  rectifying  and  non-rectifying  types.  Clearly,  because 
of  the  similarity  between  n-type  and  p-type  metal-semiconductor  con¬ 
tacts  it  has  been  a  great  temptation  to  explain  p-type  device  operation 
using  the  theory  for  n-type  devices - with,  of  course,  suitable  modi¬ 

fications.  A  particularly  difficult  extension  of  this  n-type  device  theory 
has  been  the  need  to  explain  the  reverse  current  in  p-type  Schottky  bar¬ 
riers,  when  operating  at  a  low  ambient  temperature.  At  low  temperatures 
the  hole-electron  pair  generation  is  insufficient  to  explain  the  large 
reverse  electric  current  exhibited  by  this  semiconductor  diode. 

In  reverse  biased  n-type  Schottky  barrier  a  nearly  unlimited  sup¬ 
ply  of  mobile  electrons  reside  within  the  metallic  contact.  If,  indeed, 
within  this  source  some  electrons  have  sufficient  energy  to  overcome  the 
interface  potential  barrier,  these  electrons  will  contribute  to  the  reverse 
electric  current.  Contrasting  with  this  n-type  device  concept,  no  such 
obvious  supply  of  moble  carriers  exists  for  the  p-type  Schottky  barrier. 

It  is  evident  that  the  reverse  electric  current  in  a  p-type  device  arises 
from  holes  moving  from  the  depletion  layer  to  charge  neutral  regions  of 
the  semiconductor.  An  important  question  is  the  source  of  these  holes, 
particularly  at  an  ambient  temperature  that  decreases  generation/recom¬ 
bination  to  a  negligible  value. 

It  has  been  a  generally  accepted  practice  to  assume  hole  genera¬ 
tion  within  the  metallic  contact  of  a  p-type  Schottky  barrier  and,  thereby, 
imply  a  substantial  degree  of  similarity  between  n-type  and  p-type  de¬ 
vice  operation.  If,  indeed,  holes  were  generated  in  the  metal,  the  theory 
of  p-type  diode  operation  would  be  a  simple  extension  of  n-type  device 
theory.  It  must  be  acknowledged  that  mobile  holes,  as  such,  cannot 
exist  within  a  metallic  conductor  and,  therefore,  the  theory  of  p-type 
Schottky  barrier  operation  is  in  great  difficulty. 
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There  is  little  question  that  the  reverse  electric  current  must 
arise  from  mobile  holes  within  the  space-charge  region  of  a  p-type  de¬ 
vice.  Additionally,  it  must  be  acknowledged  that  these  mobile  holes 
do  not  arise  from  hole-electron  pair  generation  alone.  If  generation/ 
recombination  mechanisms  represented  the  only  source  of  mobile  holes, 
we  could  not  explain  the  large  magnitude  of  reverse  current  at  low  de¬ 
vice  temperatures.  Instead,  new  mechanisms  of  hole  generation  must  be 
proposed;  mechanisms  whereby  mobile  holes  are  generated  without  neces¬ 
sarily  generating  conduction  band  electrons. 

The  basis  of  the  present  research  is  to  investigate  a  new  mecha¬ 
nism  of  the  foregoing  type.  Specifically,  it  is  proposed  that  hole  genera¬ 
tion  can  exist  at  the  metal-semiconductor  interface  by  valence  electron 
transitions  from  the  semiconductor  to  unoccupied  states  within  the  metal¬ 
lic  conductor.  Further,  it  is  suggested  that  through  this  mechanism  we 
can  explain  the  low  temperature  reverse  electric  current  characteristics 
exhibited  by  low-barrier  Schottky  devices  on  p-type  silicon.  Additionally 
it  is  suggested  that  this  mechanism  will  often  render  such  a  physical  struc¬ 
ture  electrically  Ohmic  at  room  temperature. 

2.0  Field  Emission  from  Semiconductor  Materials 

For  many  years  electron  emission  from  electrical  conductors  has 
been  a  topic  of  extensive  investigation.  There  are  two  basic  mechanisms 
whereby  this  emission  can  take  place:  first,  at  high  temperatures  (Thermi¬ 
onic  emission)  and,  second,  when  a  conductor  is  exposed  to  a  large 
electric  field  in  a  direction  normal  to  its  surface  (Field  emission).  At 
high  temperatures  free  electrons  within  a  conductor  acquire  sufficient 
energy  to  overcome  an  electrostatic  potential  barrier  present  at  its  surface. 
At  large  values  of  normally  directed  electric  field  the  surface  barrier  of  a 
material  can  become  exceedingly  thin  and,  thereby,  permit  tunneling  from 
the  material  to  free  space. 


Emission  experiments  using  semiconducting  material  have  brought 
renewed  interest  to  this  subject.  In  particular,  two  distinctly  different 
field-emission  characteristics  can  be  demonstrated  for  n-type  and  p-type 
silicon  and  germanium  [  1—4 3  .  Further,  some  experiments  show  that  n-type 
silicon  can  exhibit  the  field  emission  properties  of  n-type  or  p-type  sili¬ 
con.  These  two  distinctly  different  emission  characteristics  have  been 
explained  theoretically  [5]  ;  the  first  type  is  proposed  to  result  from  the 
emission  of  valence  band  electrons. 

It  is  emphasized  that  these  studies  were  directed  toward  the  topic 
of  field-emission,  with  no  reference  to  implications  associated  with  de¬ 
vice  theory.  Nonetheless,  field-emission  from  the  valence  band  repre¬ 
sents  hole  generation  without  the  generation  of  an  associated  conduction 
band  electron.  If,  indeed,  these  field-emission  workers  demonstrate  va¬ 
lence  band  emission  from  silicon,  there  is  little  reason  to  neglect  this 
mechanism  in  the  theory  of  Schottky  barrier  diode  operation. 

3 .0  Unoccupied  Energy  States  in  a  Metallic  Conductor 

In  the  present  discussion  we  focus  upon  the  density  and  distri¬ 
bution  of  unoccupied  energy  states  in  a  conductor  that  reside  between 

the  Fermi- level  E.  and  the  conduction  band  edge  E  .  It  will  be  shown 
t  c 

that  a  significant  quantity  of  these  unoccupied  states  are  far  removed 

from  the  Fermi-level,  and  there  is  no  reason  to  assume  Er  -  E,  in  Fig.l 

t  b 

represents  a  real  energy  barrier  in  a  p-type  Schottky  device.  In  partic¬ 
ular  many  silicon  valence  electrons  at  the  metal-semiconductor  interface 
that  undergo  thermal  excitation  can,  indeed  ,  enter  the  metal  at  an  energy 
substantially  below  the  traditionally  assumed  Fermi-level. 

The  distribution  of  unoccupied  energy  states  is  particularly  de¬ 
pendent  upon  the  position  of  the  Fermi-level,  relative  to  the  conduction 
band  edge.  Additionally,  this  Fermi-level  is  temperature  dependent;  we 
must  first  establish  its  magnitude,  as  a  function  of  temperature.  This 
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is  easily  accomplished  using  the  conventional  density  of  states  equa¬ 
tion  [6 1  , 

AN  =  8n\/2~ (v^n/hjVfT  AE  (1-1) 


in  combination  with  the  Fermi-Dirac  distribution 

-1 


P(E)  =  [l  +  exp  (E  -  E{)  /  kl] 


(1-2) 


Upon  combining  Eq.(l-l)  and  Eq.(l-2)  we  have  the  familiar  density  of 
occupied  energy  states  that  reside  between  E  and  (E  +  AE), 


AN 


o  r  *1“  1 

=  8n'v/2~(\Sm~/h)  VT  ll  +  exp  (E  -  Ef)  /  kT  AE  (1-3) 


clearly,  the  total  density  of  these  occupied  states  in  a  metal  is  given  by 

oo 


N  =  8fr\/2-(\/nv'hr 


V/TT  dE 


(1-4) 


1  +  exp  (E  -  E  )  /  kT 


Equation  (1-4)  has  been  used  to  calculate  the  Fermi-level  in  a  metal  as 
a  function  of  temperature.  We  first  assume  in  Eq.(l-4)  that  T  =  0  yielding 


N  =  8  v\fl~  (v/nT/h)^ 
c 


2 

3 


(1-5) 


where  N  is  the  density  of  conduction  electrons  and  E,n  is  the  Fermi 
c  fu 

energy  at  T  =  0.  Thereafter,  we  combine  Eq . (1-4)  and  Eq . (1- 5)  and  solve 
for  E{ 
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1  +  exp  (E  -  E  )  /  kT 


Figure  2  illustrates  this  Fermi-level  calculation  assuming  a  den¬ 
sity  of  states  effective  mass  of  unity.  Clearly,  calculating  the  Fermi- 
level  is  an  iterative  process  and  must  be  accomplished  through  numerical 
integration  techniques;  the  computer  program  used  for  this  calculation 
is  listed  in  APPENDIX  IA.  In  all  calculations  presented  in  this  report  the 
forementioned  computer  program  was  used  for  calculating  the  Fermi-level 
Only  a  minor  modification  of  Eq.(l-4)  is  required  to  calculate  the 
density  of  unoccupied  energy  states  that  reside  between  the  Fermi-level 
and  the  conduction  band  edge. 


N  =  8tt  \T7T  ( s/rn  /  h  ) 


■h 


viT  dE 


exp  (E  -  E^)  /  kT 


Figure  3  illustrates  this  calculation  throughout  a  wide  temperature 
range,  and  throughout  a  range  of  free  electron  density  frequently  en¬ 
countered  in  a  metallic  conductor. 

From  Fig. 3,  at  77°  K  conductors  with  a  free  electron  density 
20  23 

within  the  range  10  ^  Nc  <  10  can  be  expected  to  contain  about 

18  19 

10  to  10  unoccupied  energy  states  between  the  Fermi-level  and  the 
conduction  band  edge.  Clearly,  this  is  a  non-negligible  quantity  if 
we  wish  to  assume  in  Fig.l  represents  an  energy  barrier  that 

limits  the  reverse  electric  current  in  a  Schottky  diode. 


Density  of  Unoccupied  Energy  States  in  a  Metal  Between  the  Conduction 
Band  Edae  E  and  the  Fermi-Level.  N  =  Free  Electron  Density 


An  equally  important  aspect  of  this  concept  is  the  distribution 
of  these  unoccupied  energy  states  within  the  range  E  l  E  ^E  .  Clearly, 

C  I 

if  the  majority  of  these  unoccupied  states  are  clustered  near  the  Fermi- 

level  we  could  retain  the  concept  of  a  potential  barrier  of  E,  -  E,  in  a 

t  b 

Schottky  diode.  If,  instead,  a  large  quantity  of  states  reside  near  the 
energy  E^  in  Fig.l,  we  must  reconsider  the  concept  of  reverse  electric 
current  in  this  semiconductor  device. 

Again,  by  numerically  integrating  the  expression 


(1-8) 


We  evaluate  the  quantity  of  unoccupied  states  that  reside  throughout  the 

energy  range  E^4  E  4 E^.  Such  calculations  are  presented  in  Fig.  4 

through  Fig.  6  for  a  range  of  assumed  locations  for  E^. 

This  sequence  of  calculations  clearly  show  that  the  density  of 

unoccupied  states  are,  indeed,  clustered  near  the  Fermi-level.  For 

22  -3 

example,  material  containing  a  free  electron  density  of  10  (cm  ), 

Fig. (6),  exhibit  a  negligible  difference  in  vacant  state  density  between 

E./E,  =  0.1  and  E,  /E.  =  0.98,  at  a  temperature  of  77°  K.  Similarly,  at 

b  f  b  f  21-3 

room  temperature  material  with  a  free  electron  density  of  10  (cm  ), 

19 

Fig. 5,  exhibits  an  unoccupied  energy  state  density  of  about  1.25  x  10 

-3 

(cm  )  within  2%  of  the  range  (E^  -  E^)  from  the  Fermi-level.  None¬ 
theless,  it  is  recognized  that  in  a  Schottky  diode  the  remaining  unoc¬ 
cupied  energy  states  could  produce  a  significant  magnitude  of  reverse 
current . 

A  further  amplification  of  this  topic  is  obtained  by  numerically  in¬ 


tegrating  Eq.(l-4)  throughout  the  range  E  4  E  4  E  , 

c  b 


EB/EF 

8.18 


an  Energy  E,  where  E„  <  E,  <  E, 


(1-9) 


This  calculation  yields  the  density  of  unoccupied  states  that  reside  in 

the  lower  energy  regions  of  the  conduction  band. 

Figure  7  through  Fig. 10  illustrate  this  calculation  throughout  a 

20  23 

range  of  free  electron  density  of  10  4  Nc  410  .  These  calculations 

clearly  show  that  material  containing  a  larger  free  electron  density  have 
the  smallest  density  of  unoccupied  states  in  the  lower  regions  of  the 

conduction  band.  For  example,  at  71°  K  a  free  electron  density  of 

20  -3  14  -3 

10  (cm  ),  Fig. 7,  exhibits  about  10  (cm  )  unoccupied  states  in  the 

lower  50%  between  E.  and  E  .  At  this  same  temperature  material  con- 
*  c  22-3 

taining  a  free  electron  density  of  10  (cm  )  contain  approximately 
8  -3 

10  (cm  )  unoccupied  states  in  the  lower  92%  of  the  energy  between 

E,  and  E  . 
f  c 

The  reason  for  this  foregoing  situation  is  attributable  to  the  in¬ 
crease  of  (E,  -  E  )  with  free  electron  density.  Fig. 2.  An  increase  of 
f  c 

£  -  E  produces  a  decrease  of  probability  that  the  lower  energy  electrons 
X  c 

can  be  excited  to  an  energy  level  above  E^. 

In  the  following  sections  of  this  report  the  computer  programs  used 
for  calculations  thus  far  illustrated  are  an  integral  part  of  all  calculations 
directed  toward  the  operation  of  Schottky  barrier  devices . 
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Conduction  Be 


Edge  and  an  Energy  E.  where  V .  <  E^  <  Ef 


Fig. 9  Density  of  Unoccupied  Energy  States  Between  the  Conduction  Band  Edge 


4.0  The  Transition  Energy  in  a  Schottkv  Barrier 


In  the  conduction  band  of  a  metal  the  density  of  unoccupied 
energy  states  between  £  and  (E  +  AE)  is  given  by 


-i 

AN  =  8 rrV2~ ( /  h)3  n/eT  j^l  +  exp  (Ef  -  E)  /  kT^  AE  . 


(1-10) 


If,  indeed,  semiconductor  valence  band  electrons  are  to  enter  this  con¬ 
duction  band  there  must  be  a  place  for  them  to  go.  If  the  number  of 
electrons  attempting  to  enter  this  conduction  band  at  an  energy  E  ex¬ 
ceed  the  density  of  unoccupied  states,  many  could  return  and  not  con¬ 
tribute  to  the  electric  current.  If,  instead,  the  density  of  unoccupied 
states  at  an  energy  E  exceeds  the  source  of  supply,  all  valence  band 
electrons  attempting  to  enter  the  metal  conduction  band  will  find  a  place 
to  go,  and  will  contribute  to  the  barrier  electric  current. 

Next  we  consider  the  density  of  semiconductor  valence  band 
electrons  that  are  attempting  to  enter  this  metallic  conduction  band . 

From  Fig.l,  all  valence  band  electrons  at  an  energy  in  excess  of  E^  must 
have  undergone  thermal  excitation.  Because  all  energy  eigenvalues  of 
the  wave  equation  become  imaginary  in  the  forbidden  band,  we  assume  a 
Boltzmann  distribution  adequately  describes  the  density  of  electrons  in 
this  region  of  the  structure  [  6  ]  . 

1/2 

AN  =  N  exp  (-  \  Av  .  (1-11) 

v  \  2  k  T/  \  kT  /  X 

Here  we  describe  the  density  of  valence  band  electrons  with  an  x-axis 
directed  velocity  between  v^  and  (v  +  Av^)  that  are  attempting  to  enter 
the  metallic  conduction  band . 
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Clearly,  if  the  density  of  unoccupied  states,  Eq.  1-10, ex¬ 
ceeds  the  density  of  excited  valence  electrons,  Eq.1-11,  all  can 
enter  the  conduction  band.  If,  instead,  the  available  electrons 
with  a  velocity  v^  exceed  the  number  of  available  unoccupied 
energy  states, many  will  return  to  the  valence  band.  For  the  pur¬ 
pose  of  this  analysis  we  define  a  transition  energy  E  .  The  trans¬ 
ition  energy  is  that  energy  bounding  these  two  forementioned 
energy  regions .  At  energies  above  E  the  rate  at  which  electrons 
enter  the  conduction  band  is  determined  by  the  Maxwellian  distri¬ 
bution  of  valence  electrons,  Eq.1-11.  At  energies  below  E^  the 
rate  at  which  electrons  enter  the  conduction  band  is  determined 
by  the  availability  of  unoccupied  energy  states. 

Figure  11  through  Figure  14  illustrates  the  calculated  value 

of  E^,  described  in  terms  of  a  barrier  energy  =  (E^  -  E^) . 

Because  this  transition  barrier  energy  is  not  unlike  the  classical 

Schottky  interface  energy  barrier  W  where  W  =  E ,  -  E,  ,  we 

c  c  t  b 

have  related  the  two  in  these  illustrative  calculations.  Classical 
theory  is  based  upon  an  assumption  that  excited  valence  electrons 
must  reach  an  energy  E^  to  enter  the  conduction  band;  this  is  equiv¬ 
alent  to  the  general  concept  that  a  hole  at  E^.  moves  to  E^.  Con¬ 
trasting  with  this  classical  concept  the  transition  barrier  Wc  repre¬ 
sents  the  energy  to  which  valence  electrons  must  be  excited  for  all 
to  enter  the  conduction  band.  Below  W  their  entry  into  the  con¬ 
duction  band  is  governed  by  the  availability  of  unoccupied  energy 
states. 

From  these  illustrations,  Fig.  11  through  Fig.  14,  it  is  clear  that 
the  transition  energy  W  is  significantly  lower  than  the  classical  barrier  W  , 
for  all  values  of  free  electron  density.  It  is  through  calculations  of  the 
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Fig.  11  Calculated  Transition  Barrier  W  in  a  p-type  Schottky  Jevice 


Calculated  Transition  Barrier  W  in  a  p-type  Schottky  Device 


Calculated  Transition  Barrier  W  in  a  p-type  Schottky  Device 
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barrier  electric  current  above  and  below  the  transition  energy  E^  we  shall 
establish  the  reverse  electric  current  in  a  p-type  Schottky  barrier. 

5.0  The  Electric  Current  in  a  Schottkv  Barrier - Classical  Concept 

The  classical  concept  of  p-type  Schottky  diode  operation  is  that 
the  total  electric  current  can  be  represented  by  an  expression  of  the  form 

J  =  Jq  £l  -  exp  {-  q  Vq  /k  T)j  ,  (1-12) 

where  V  is  an  applied  voltage  and  J  is  the  saturation  current, 
a  o 

To  date,  a  rigorous  calculation  for  J  has  not  been  found  in  the 

o 

technical  literature.  It  is  an  accepted  practice  to  attribute  Jq  to  two  dis¬ 
tinctly  different  mechanisms:  first,  hole-electron  generation  in  the  p-type 
semiconductor  material  and,  second,  hole  transitions  from  the  metal  to  the 
semiconductor  material.  The  first  mechanism  can  be  acknowledged  as  well 
founded  but,  instead,  there  remains  no  real  quantitative  explanation  for 
the  saturation  current  at  low  temperatures;  at  temperatures  that  render  hole- 
electron  generation  insignificant. 

In  order  to  adequately  evaluate  our  concept  of  p-type  device  opera¬ 
tion  it  is  necessary  to  recast  this  classical  theory  into  a  form  suitable 
for  comparison.  For  this  reason  we  abandon  the  notion  of  hole  generation 
in  the  metal  and  consider  this  traditional  theory  in  terms  of  valence  elec¬ 
tron  excitation  to  the  Fermi-level.  Clearly,  this  is  consistant  with  the 
concept  of  hole  transitions  from  the  Fermi-level  to  the  valence  band. 

As  in  Section  4.0  of  this  chapter,  we  assume  a  Maxwellian  ve¬ 
locity  distribution  for  valence  electrons  given  by 

AN  =  N  ( — — — ^  exp  /-  -i—  invVkl)  Av  (1-13) 

v  y  2  ffk  T  /  V  2  x  /  x 
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Further,  we  assume  an  electric  current  produced  by  this  velocity  dis¬ 
tribution  given  by 


where 


v 


f 


Vf  =  (2  E  /  k  T)  . 


(1-15) 


Equation  (1-14)  represents  an  expression  for  the  saturation  cur¬ 
rent  J  that  retains  the  traditional  concept  of  hole  transitions  from  the 
o 

Fermi-level  to  the  conduction  band.  Nonetheless,  this  expression  has 

been  developed  in  a  manner  that  we  believe  is  quantitatively  reasonable, 

and  can  be  used  for  comparison  with  new  concepts  of  device  operation 

arising  from  the  present  research. 

In  the  foregoing  equations  represents  the  density  of  valence 

electrons  that  are  at  an  energy  in  excess  of  the  valence  band  edge  at  the 

metal-semiconductor  interface,  E^.  This  is  easily  calculated  if  we  first 

determine  the  location  of  E  ,  relative  to  the  lower  edge  of  the  valence  band. 

v 

The  density  of  occupied  energy  states  in  a  semiconductor  is  given  by 

o 


vTT  dE 


1  -  exp  (E  -  E  )  k  T 


.  (1-16) 
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If  we  assume  a  temperature  of  zero  degrees  K,  the  Fermi-level  is  es¬ 
sentially  at  Ev<  and  we  have 


N 


t 


8n  \fl 


(1-17) 


where  N  is  the  total  number  of  valence  electrons,  and  E,„  is  the  Fermi- 

t  fo 

level  (at  0°  K)  as  measured  with  respect  to  the  lower  boundary  of  the 

valence  band.  Because  at  the  interface  E.n  and  E,  coincide,  we  can  now 

f  0  b 

state  that  the  density  of  unoccupied  states  in  the  valence  band  is 


N  = 
v 


\Ze~  dE. 

1  +  exp  (EJ  +  E)  /  k  T 


(1-18) 


where  E'  is  the  interface  Fermi-level  when  measured  with  respect  to  the 
lower  edge  of  the  semiconductor  conduction  band. 


e;  -  Eto + <Ef  -  y 


(1-19) 


Because  valence  electrons  not  found  in  the  valence  band  must 

be  at  an  energy  in  excess  of  E,  ,  we  know  the  number  of  electrons  in- 

fa 

volved  in  the  Maxwellian  distribution  and,  hence,  for  Eqs.(l-12)  and  (1-13). 

6.0  Reverse  Current  From  Transitions  Above  E^ 

As  outlined  in  Section  4.0,  electron  excitation  from  the  valence  band 
E^,  Fig.l,  to  an  energy  above  E  (the  transition  energy)  will  assure  a  trans¬ 
ition  into  the  metal.  Hence,  like  the  traditional  theory,  we  can  immediately 
establish  this  component  of  the  diode  saturation  current  through  a  minor 
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modification  of  Eq.1-14. 


=  q  I  v^dn  = 


i. 

v  (t^t)  axp  (-EtAT) 


(1-20) 


where 


v  =  (2  Et  /  m) 


(1-21) 


This  expression  has  been  used  in  all  of  the  following  calculations  as 
one  component  of  the  total  saturation  current  J  ,  and  is  determined 
using  Eq.1-18. 


7 . 0  Reverse  Current  From  Transitions  Below  Et 

Below  the  transition  energy  Et  the  density  of  unoccupied  energy 
states  becomes  less  than  the  number  of  electrons  available  to  occupy  these 
states.  If  we  recast  the  density  of  unoccupied  states  into  an  expression 
for  electron  velocity,  Eq.1-10,  we  obtain 


AN  =  Nc  — y  (Ef[)  /  m)  3^2  v2  £l  +  exp  (E{-  y  mv2)J  dv  (1-5 


Thus,  the  electric  current  produced  by  transitions  to  these  energy  states 
is  given  by 


,2  Nc  f  (Ef(/m)-3/2 


3 

_ v  d  v _ 

1  +  exp  (Ef  “  —  v) 


(1-23) 
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where 


1 

vt  =  (2  Et  /  m)"2" 


(1-24) 


Rewriting  Eq.1-19  into  energy  we  obtain 


E  dE 

1  +  exp  (E  -  E)  /  k  T 


(1-25) 


Equation  (1-21)  is  the  expression  we  use  for  calculating  the  electric 
current  produced  by  electron  transitions  below  the  transition  energy  E^. 


8-0  Computations  of  Reverse  Electric  Current 


Few,  if  any,  of  the  foregoing  equations  can  be  evaluated  on  an 
analytical  basis;  integration  of  the  Fermi-Dirac  distribution  equation  must 
be  done  numerically,  on  a  computer.  For  this  reason,  our  initial  evalua¬ 
tion  of  this  new  concept  of  Schottky  barrier  operation  has  been  accomplished 
through  a  sequence  of  electric  current  calculations.  Equation  (1-20)  and 
Eq.(l-25)  have  been  used  for  calculating  the  current  predicted  by  these 
new  theoretical  concepts,  and  such  calculations  are  compared  with  those 
using  Eq.1-14  as  representative  of  the  classical  theory  for  Schottky  diode 
operation . 

Figure  15  through  18  illustrate  the  total  reverse  current  predicted 
by  Eq.(l-20)  and  Eq.(l-21)  for  J  and  J  respectively,  throughout  a  wide 
assumed  range  for  the  free  electron  density  in  the  metal  conduction  band. 

A  comparison  is  first  undertaken  between  these  calculations  Fig.  15  through 
Fig.  18,  and  similar  calculations  based  upon  classical  theory.  Fig. 19  and 
Fig. 20.  Because  W  is  given  in  these  calculations  (W  =E,-E  ),  classical 
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theory  is  nearly  independent  of  the  free  electron  density  in  the  metal  and, 

hence,  only  a  single  illustration  is  needed  for  J  .  For  convenience,  J 

c  c 

is  presented  in  two  illustrations.  Fig.  19  and  Fig.  20,  which  have  different 

assumed  ranges  for  W  . 

c 

We  encounter  two  different  situations  in  these  calculations  with 

relation  to  N  ,  the  conduction  band  electron  density:  first,  when  N  is 

C20  23  C 

small  (N  slO  )  and,  second,  when  N  is  large  (N  arlO  ).  When  N 
c  c  c  c 

is  small  the  metallic  conduction  band  closely  resembles  that  of  a  highly 
degenerate  n-type  semiconductor;  the  Fermi-level  is  only  a  small  dis¬ 
tance  above  the  conduction  band  edge,  E^.  As  a  consequence,  many  of 
the  assumed  values  for  barrier  energy  W^  places  the  semiconductor  valence 

band  edge,  E,  ,  below  the  conduction  band  edge,  E  .  In  the  second  case 

23  b  c 

(N  SlO  )  we  have  a  metallic  conductor  in  which  the  Fermi-level  is  well 
c 

above  the  conduction  band  edge  and,  therefore,  any  reasonable  assump¬ 
tion  for  W  (W  4  2: 1.0  ev)  places  the  semiconductor  valence  band  edge 
c  c 

at  the  interface,  E,  ,  well  within  the  metal  conduction  band  (E  4  E,  4  EJ. 

b  20  c  d  { 

Calculations  show  that  by  assuming  N  =  10  ,  this  new  theory 

c 

of  operation  yields  a  reverse  current,  J  ,  that  is  always  substantially 
larger  than  that  given  by  classical  theory,  J  .  For  purposes  of  compari¬ 
son,  some  illustrative  calculations  yield  at  T  =  60°  K: 


W  (ev) 
c 

2 

Jt  (amps/cm  ) 

2 

Jc(amps/cm  ) 

0.5 

10-53 

-60 

10 

0.3 

10-23 

io-29 

0.25 

10-15 

10-22 

0.2 

10-8 

io-14 

Because  in  the  foregoing  calculations  lies  below  the  conduction  band 

edge  no  contribution  to  J  arises  from  direct  electron  transitions  from  the 

valence  band  to  the  conduction  band;  all  transitions  are  due  to  thermally 

20  -3 

excited  valence  electrons  with  an  energy  in  excess  of  E,  .  For  N  =  10  (cm  ), 

b  c 
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only  at  an  assumed  barrier  of  .05  ev  (W  =  0.05)  does  this  direct  transi- 

c 

tion  mechanism  produce  an  electric  current,  J  ,  that  contributes  signifi¬ 
cantly  to  J  ,  Fig. 25.  For  this  low  barrier  situation  (W  =0.05  ev)  we 

L  C 

calculate: 
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Contrasting  with  the  foregoing  small  free  electron  situation,  by 
23  -3 

assuming  N  =10  (cm  )  the  Fermi-level  is  moved  well  away  from  the 
C  20 

conduction  band  edge.  As  with  Nc  =  10  we  again  obtain  a  small  elec¬ 
tric  current  component  due  to  direct  transitions ,  but  for  an  entirely  dif- 

23 

femt  reason.  Setting  N  =10  places  the  interface  valence  band  energy, 

c 

E,  ,  well  within  the  metallic  conduction  band  edge  which,  under  the  proper 
b 

circumstances,  could  produce  a  large  direct  transition  current.  When 
23 

=  10  we  find  that  unoccupied  conduction  band  energy  states  are 

clustered  within  a  small  energy  range  from  E^.  As  a  consequence,  when 

E  is  sufficiently  large  we  obtain  a  negligible  amount  of  current  due  to 
c 

direct  electron  transition  from  the  valence  band  to  the  conduction  band. 

For  the  purposes  of  illustration,  the  following  represent  the  calcu¬ 
lated  electric  current  at  T  =  83°  k: 


r  (ev) 
c 
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Saturation  Current  Component  Due  to  Transition  of  Electrons 
with  an  Energy  Greater  Than  E^ 


Fig .  21  Saturation  Current  Component  Due  to  Transition  of  Electrons 
with  an  Energy  Greater  Than  E 


Saturation  Current  Component  Due  to  Electrons  with  an 
Energy  Greater  Than  E^ 


Saturation  Current  Component  Due  to  Transitions  by  Electrons 
with  an  Energy  Less  Than  E^ 


where  current  is  in  amps/cm  .  In  all  of  these  calculations  the  direct 

-47  2 

transition  current  remained  relatively  unchanged  at  10  amps/cm  ; 
this  situation  probably  arises  from  underflow  protection  within  the  com¬ 


puter  program.  This  sequence  of  calculations  indicate  that  very  few  unoc- 

23 

cupied  states  reside  in  lower  regions  fo  the  conduction  band  when  =  10 
Nonetheless,  a  sufficient  number  exist  below  so  that  J  is  significantly 
larger  than  J  . 

c  20  23 

An  intermediate  range  is  found  between  N  =10  and  N  =10 

c  c 

where  both  J.  and  J  can  contribute  in  a  significant  fashion  to  the  total 

&  20  23 

electric  current.  In  this  intermediate  region,  where  D  <  N  <  10  , 

c 

the  unoccupied  energy  states  can  reside  at  energy  levels  far  below  the 

Fermi-level,  E^.  As  a  consequence  of  this  situation  the  transition  energy 

E  can  reside  well  below  Er  and,  therefore,  W  <  W  ,  which  assures  J  <  J, . 
t  t  t  c  cl 

Additionally,  a  large  number  of  unoccupied  energy  states  can  reside  below 
the  transition  barrier  W^,  and  this  situation  implies  a  direct  transition  current 
J2  that  can  be  significant. 

We  first  consider  the  electric  current  calculated  for  an  assumed 

21  -3 

free  electron  density  of  10  (cm  ).  At  66°  K  our  calculations  yield  the 
following  results: 


W  (ev) 
c 

Jt 

J1 

J2 

Jc 
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00 

3 

where  all  current  is  in  amps/cm  .  In  this  situation  two  different  (although 

related)  mechanisms  contribute  to  the  total  electric  current  J  .  Transitions 

of  thermally  excited  electrons  to  unoccupied  energy  states  produce  an 

electric  current  component  that  is  significantly  greater  than  predicted 

from  classical  concepts,  J  .  This  situation  implies  that  the  transition  to 

c 

unoccupied  energy  states  below  E^  (E^  4  E  4E^)  occurs  at  a  greater  rate 
than  transition  to  unoccupied  states  above  E^  (E^4  E  4»).  These  calcula¬ 
tions  also  show  that  transitions  with  no  apparent  thermal  excitation 
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Saturation  Current  Component  Due  to  Transitions  by 
Electrons  with  an  Energy  Less  Than  E 


(E^  E  4  E  )  take  place  at  a  rate  that  is  significantly  greater  than  the 
transition  of  thermally  excited  valence  electrons.  Hence,  in  this  com¬ 
putational  illustration  both  and  are  significant  in  magnitude,  and 
in  combination  they  produce  a  greater  electric  current  than  predicted 

using  classical  concepts  of  diode  operation. 

22  -3 

By  assuming  N  =10  {cm  )  we  obtain  a  Fermi-level  that  is  well 
c 

above  the  conduction  band  edge,  yet  the  distribution  of  unoccupied  states 
extends  for  a  significant  energy  range  below  this  Fermi-level.  As  a  con¬ 
sequence,  direct  transition  current  can  be  very  large.  Fig.  26,  in  con¬ 
junction  with  a  large  current  due  to  transitions  above  the  transition 
energy  E^,  Fig.  23.  In  combination,  the  total  junction  current,  Fig.  17, 
is  much  larger  than  expected  from  conventional  concepts  of  device  opera¬ 
tion,  Fig.  19  and  Fig. 20. 

It  should  be  mentioned  that  the  apparent  change  of  slope  in  Fig. 17 
appears  to  be  a  real  mechanism,  and  not  a  problem  of  the  computer  pro¬ 
gram.  At  high  temperature  the  total  current  is  dominated  by  direct  valence 
electron  transitions  to  the  conduction  band  and,  thereby,  determining  the 
total  current  J^.  A  decrease  of  temperature  produces  a  rapid  fall  off  of  this 
direct  transition  current.  Fig. 27,  and  J  is  determined  by  transitions  above 
the  transition  energy  E^.  Although  it  is  not  shown  in  these  calculations, 
at  a  lower  temperature  many  of  these  curves  in  Fig.  17  exhibit  another 
change  of  slope;  this  change  is  attributed  to  a  rapid  decrease  of  N^;  the 

valence  electrons  residing  above  the  valence  band  edge,  E,  ,  at  the  metal- 

b 

semiconductor  interface. 

9 . 0  Reverse  Current  in  a  p-tvoe  Schottkv  Barrier  at  77  °K 

In  reverse  bias  situations  it  is  presumed  that  holes  are  generated  at 
the  metal- semiconductor  interface  through  the  foregoing  mechanisms.  After 
generation  these  holes  move  to  charge  neutral  regions  of  the  p-type 
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.26  Saturation  Current  Component  Due  to  Transitions  by 
Electrons  with  an  Energy  Less  Than  E 


material,  under  the  influence  of  a  space-charge  electric  field;  there¬ 
by  a  reverse  electric  current  is  produced  in  this  semiconductor  device. 
Experiment  shows  that  the  magnitude  of  this  reverse  current  increases 
with  applied  reverse  voltage,  and  this  mechanism  has  not  yet  been  in¬ 
cluded  in  the  present  theory. 

From  conventional  concepts  of  device  operation  it  can  be  shown 
[  6  ]  that  a  lowering  of  the  potential  barrier  results  from  the  attractive 
forces  between  an  interface  generated  hole,  and  an  image  charge  at  the 
metal.  This  mechanism  is  frequently  called  the  Schottky  effect,  and 
through  this  mechanism  E^,  Fig.l,  will  be  increased  by  an  amount 


(1-26) 

where  f q  represents  the  space-charge  layer  electric  field  at  the  metal  - 
semiconductor  interface.  Using  the  depletion  approximation  for  deter¬ 
mining  this  electric  field,  Eq.1-26  has  been  added  to  assumed  values  of 
E^  yielding  a  conventional  interface  energy  barrier  of 

wc * Ef -  "b  - -r  (1'27) 

We  next  apply  our  new  theoretical  concepts  to  a  reverse  biased 
Pt-Si  diode  operating  at  77°  K.  Additionally,  we  compare  the  results  of 
such  calculations  with  similar  calculations  using  conventional  theoretical 
concepts.  It  is  emphasized  that  little  is  known  about  the  Pt-Si  in  this 
semiconductor  device  and,  therefore,  unknown  parameters  must  be  in¬ 
troduced  by  guess,  rather  than  fact.  Specifically,  for  this  calculation 
we  require  the  free  electron  density  in  the  Pt-Si  region,  and  the  density 
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of  states  effective  mass  of  these  free  electrons.  For  the  latter  we 
assume  ( m*/m )  is  equal  to  unity,  which  is  indeed  incorrect,  but  this 
parameter  is  unknown. 

For  the  free  electron  density  we  have  adjusted  its  magnitude  to 

yield  a  reverse  electric  current  typical  of  that  obtained  in  experimental 

devices  operating  at  77°  K.  Our  best  guess  is  that  the  Pt-Si  region  con- 

20  -3 

tains  a  free  electron  density  of  approximately  10  (cm  ),  which  im¬ 
plies  this  region  is  more  accurately  classified  as  a  degenerate  semicon¬ 
ductor,  rather  than  a  metal.  Figure  27  illustrates  the  calculated  reverse 

electric  current  vs.  applied  voltage,  assuming  a  free  electron  density 
20  -3  20  -3 

between  10  (cm  )  and  1.5x10  (cm  ),  assuming  an  energy  barrier  W 

c 

of  0.27  (ev).  Similarly,  using  conventional  theoretical  concepts  the 
calculated  reverse  current  is  illustrated  in  Fig.  28;  this  calculated  cur¬ 
rent  is  independent  of  the  Pt-Si  characteristics . 

These  calculations  suggest  that  the  reverse  current  of  a  p-type 
Pt-Si  diode  will  be  many  orders  of  magnitude  larger  than  predicted  from 
traditional  theoretical  concepts.  Further, implied  by  these  calculations 
is  that  without  good  control  over  the  fabrication  processes  poor  reproduct- 
ability  might  be  expected  in  Pt-Si  devices.  For  example,  our  calculated 

reverse  current,  at  a  given  biasing  voltage,  changes  by  nearly  two  orders 

20  -3  20 

of  magnitude  for  a  free  electron  density  between  10  (cm  )  and  1.5x10 
(cm  3). 

It  can  be  shown  that  an  assumed  barrier  W  of  0.27  ev  locates  the 

c 

valence  band  edge  at  the  interface,  E^  in  Fig.l,  below  the  conduction  band 

edge,  E  .  For  this  reason  the  mechanism  of  direct  electron  transition  can- 
c 

not  take  place  between  the  valence  band  and  the  Pt-Si  region  of  this  de¬ 
vice.  Nonetheless,  calculations  show  that  the  transition  energy  lies 

close  to  the  conduction  band  edge  E  yielding  a  true  energy  barrier  of  about 

c 

E  -  E,  ,  rather  than  E  -  E.  =  0.27  (ev),  as  suggested  by  conventional 
c  b  ID 

theory.  It  is  suggested  that  this  significant  difference  between  the  assumed 
energy  barrier  and  the  true  energy  barrier  produces  the  unexpected  large 
reverse  electric  current  in  a  Pt-Si  diode. 
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Calculated  Volt-Ampere  Characteristics  Based  Upon  Results  from  this  Research 
(Assumed  Barrier  W  =  0.27  ev,  T  =  77°K) 


Calculated  Volt-Ampere  Characteristics  of  a  Schottky  Barrier  Based  upon  Con 
ventional  Device  Concepts  (Assumed  Barrier  W  =  0.27  ev,  T  =  77°  K) 


Chapter  II 


Potential  Calculations  for  a  Schottkv  Barrier  Diode 


1.0  Introduction 


Most  mathematical  models  of  semiconductor  device  operation 
are  based  upon  the  concepts  and  mathematical  equations  outlined  by 
W.  Van  RoosbroekC  7  ]in  1950.  In  one  spatial  dimension,  the  steady- 
state  form  of  these  equations  are  given  by: 


a) 


b) 


c) 
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d) 


e) 


d  J 
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=  C. 


=  c. 


f) 


JT  =  Jn  +  Jn 

T  p  n 
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Equation  (2-la)  is  Poisson's  equation,  and  it  relates  the  diver¬ 
gence  of  the  electric  field  to  the  electrostatic  charge  distribution  aris¬ 
ing  from  both  mobile  charge  carriers  (holes  p  and  electrons  n) ,  and  im¬ 
mobile  ionized  impurity  atoms  within  the  semiconductor  lattice. 

Equations  (2-lb)  and  (2-lc)  give  the  electric  current  density  aris¬ 
ing  from  the  transport  of  mobile  holes  and  electrons.  These  equations 

express  the  dependency  of  each  electric  current  component,  J  and  J  , 

P  ti 

upon  the  concentration  gradient  of  holes  and  electrons,  the  mobility 
of  these  charge  carriers,  and  the  electric  field  within  the  semicon¬ 
ductor  material. 

Equations  (2-ld)  and  (2-le)  are  the  continuity  equations  for  holes 
and  electrons  in  a  semiconductor,  and  these  equations  contain  unspeci¬ 
fied  rates  of  carrier  generation  and  recombination. 

Equation  ( 2— If )  states  that  the  total  electric  current  density,  J^, 

is  the  sum  of  the  electric  current  due  to  holes,  J  ,  and  electrons,  J  . 

P  ^ 

Throughout  this  analysis  we  assume  the  applicability  of  Boltzmann 
statistics,  and  the  hole  and  electron  densities  are  given  by 


a) 

b) 


n  =  n,  exp 


p  =  n.  exp 
1 


-  [q  (  V  -  <t>n  )  /  kTj 

-  |q  (  <*>p  -  V  )  /  kT  j 


(2-2) 


where  <t>  and  <£  are  the  quasi-Fermi  potentials  for  electrons  a^d  holes, 
n  p 

respectively . 

The  carrier  generation/recombination  term  in  Eqs.(2-ld)  and  (2-le) 
is  assumed  composed  of  three  terms:  first,  carrier  generation/recombination 
through  defects,  second  by  impact  and/or  avalanche  ionization  and,  third, 
hole  generation  through  valence  electron  transitions  from  the  valence  band 
of  p-type  silicon  to  the  metal. 


i 
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Hole/electron  generation  and  recombination  through  defects  are 
represented  by  a  Hall-Schockley-Read  (HSR)  single  level  model  [8,9] 
which  characterizes  defects  with  neutral  and  single  charge  states. 
Carrier  generation/recombination  through  a  single  defect  level  center  Do] 
is  given  by: 


R  =  C 


2 

pn  -  n. 
1 


p -  "  J. 

n+p+2n;  Cosh  q ( <f> ^  -  0.)  /  k  T 


(2-3) 


where  <f>  represents  the  trap  potential  and  ;  represents  the  intrinsic 
Fermi  potential  in  silicon. 

Calculations  suggest  that  at  the  temperatures  considered  here 
(77°  K-110”  K)  terms  in  Eqs.2-1  involving  electrons  in  p-type  silicon  will 
have  a  negligible  influence  upon  the  results  attained .  For  this  reason 
such  terms  have  been  neglected  and,  thereby,  these  equations  are  re¬ 
duced  to  the  form 


a) 


d  i/»  q 

- 5 - - 

dx  K« 
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-  p) 


b) 


dp 

dx 


- 


d^ 

dx 


c)  djP  =  dJn  =  c 

dx  dx 


(2-4) 


2.0  Electric  Current  Continuity 

Many  investigations  of  diode  operation  neglect  the  influence  of 
electric  current  upon  the  electrostatic  potential  distribution.  Such  a 
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simplification  is  achieved  by  maintaining  constant  the  electron  and 
hole  quasi-Fermi  potentials.  In  the  presence  of  a  large  ionized  impur¬ 
ity  atom  density  little  error  results  from  this  simplification,  particularly 
when  the  reverse  current  is  relatively  small.  Nonetheless,  velocity 
limited  carrier  transport  in  a  reverse  biased  diode  does  produce  a  sig¬ 
nificant  amount  of  mobile  carrier  accumulation  and,  thereby,  can  modify 
the  space-charge  layer  potential  distribution. 

For  purposes  of  this  analysis  we  attain  electric  current  continuity 
by  relating  the  divergence  of  electric  current  to  the  density  of  mobile  car¬ 
riers  at  a  given  location: 


=  R(p). 


(2-5) 


Both  drift  and  diffusion  mechanisms  are  assumed  to  contribute  to  the  total 
electric  current  in  Eq.2-5  thus: 


dp 

l  =  -  qo 

p  p  dx 


~  q  **p  P 


(2-6) 


3 .0  Temperature  Dependent  Parameters 

At  low  temperatures  (77°K-110°K)  many  physical  and  electrical 
properties  of  silicon  differ  from  that  realizes  at  room  temperature.  Un¬ 
fortunately,  substantial  work  remains  to  accurately  characterize  silicon 
at  low  temperature  but,  nonetheless,  it  is  believed  that  many  important 
low  temperature  material  parameters  are  available  in  the  technical  literature. 


The  energy  gap  Eg  in  silicon  can  be  shown  slightly  tempera¬ 
ture  dependent Ql  ] ,  and  approximated  by  the  relation: 


Eg  =  1.21  -  7 . 242xl0~3  (T/300)-3 .664x10“ 2(T/300)2 


(2-7) 
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This  expression  represents  an  empirical  fit  to  Fig. 2-10 
of  Ref  .11 . 

b)  Intrinsic  Fermi-Level 

The  intrinsic  Fermi-level  Ei  is  usually  assumed  located 
at  the  mid-gap  of  silicon.  Nonetheless,  in  more  rigorous 
calculations  the  intrinsic  Fermi-level  deviates  slightly  from 
this  midgap  position: 

E.-Ev  =  Eg/2  +  (k  T/q)  In  (N yu  ) 

=  Eg/2  -  1.206x10_2(T/300),  (2-8) 

where  Nv  and  Nc  represent  the  effective  density  of  states 
at  the  valence  and  conduction  band  edges,  respectively. 

In  Eq.2-8  we  assume  (N ^/N^)  =  1.02/2.8,  as  stated  in 
Ref. 11,  p.359 . 

c)  Intrinsic  Carrier  Density 

The  intrinsic  carrier  density  in  silicon  is  exceedingly 
dependent  upon  the  temperature,  as  shown  in  Ref .12,  p.27: 

n2  =  4.9xl015(M*  M  /M2)  3//4  T3^2  exp(-E  /2  k  T)  (2-9) 

l  he  g 

•jf  ^ 

where  and  Mg  are  the  density  of  states  effective  mass 
for  holes  and  electrons,  respectively,  and  M  is  the  free 
electron  mass.  This  expression  is  approximated  by 

19  3/2 

n  -  3.925x10  (T/300)  '  exp  (-  E  /2  k  T)  (2-10) 

In  Eq.2-I0  the  leading  constants  have  been  adjusted  so  that 
10  3 

ni  equals  1.5x10  /cm  when  T  =  300°  K. 
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d)  Bulk  Hole  Density 

In  charge-neutral  regions  of  p-type  silicon  the  hole  den¬ 
sity  is  given  by 


p  =  n.  exp  (E.  -  Ef)  /  k  T  (2-11) 

where  n^  represents  the  intrinsic  carrier  density,  E^  is  the 
intrinsic  Fermi-level,  and  is  the  Fermi-level.  Clearly, 
the  hole  density  in  charge  neutral  p-type  silicon  is  exceed¬ 
ingly  dependent  upon  the  temperature . 

e)  Ionized  Acceptor  Density 

The  ionized  acceptor  density  can  be  calculated  using  an 
expression  of  the  form[l3]: 

na 

N  =  - - -  (2-12) 

1  +  ^4+2  exp  {-qA/k  T)J  exp  [-q(EA  -  Ep)  /  k  tJ 

where  A  =  .044  ev;  this  term  accounts  for  the  spin  orbit  split¬ 
ting  at  the  valence  band  maximum  in  silicon . 

f)  Fermi- Level 

Calculation  of  the  Fermi-level  E^,  at  a  given  temperature 
represents  an  iterative  process.  Specifically,  at  any  given 
temperature  T  the  Fermi-level  must  be  such  that  charge  neu¬ 
trality  exists  in  the  material  under  consideration: 

-N  +  p-n  =  0.  (2-13) 

In  this  calculation  we  ignore  the  electron  density  since  at 
low  temperatures  this  term  has  negligible  influence  upon 
Eq.2-13.  Substituting  Eq.2-10  into  Eq.2-11  yields  the  hole 
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density  in  terms  of  the  Fermi-level  E^.  Thereafter,  we  sub¬ 
stitute  this  expression  and  Eq.2-12  into  Eq.2-13.  We  have 
the  electrostatic  charge  in  intrinsic  silicon,  as  a  function 
of  the  Fermi-level  E^.  Thereafter  the  Fermi-level  can  be 
established  using  iterative  methods  to  produce  an  electro¬ 
static  charge  of  zero. 

Although  the  low-temperature  mobility  for  holes  is  available  in  the 
technical  literature  [l3],  to  date,  this  has  not  been  implemented  into  this 
mathematical  model. 


4 . 0  Finite-Difference  Implementation 

Using  traditional  finite-difference  techniques,  this  problem  is 
solved  on  a  one-dimensional  grid  containing  a  specified  number  of  points 
(MAX  =  100).  These  points  are  assumed  equally  spaced,  and  located  be¬ 
tween  the  ohmic  contact  to  p-type  silicon,  and  the  metal- semiconductor 
interface. 

In  one  spatial  dimension,  we  designate  a  unit-cell  of  this  one¬ 
dimensional  array  as  composed  of  three  locations:  left,  right,  and  center. 
Thereafter,  applying  central  difference  techniques  in  conjunction  with  the 
Taylor  series,  we  can  express  the  second  derivative  of  electrostatic  po¬ 
tential  at  this  center  location  by  the  relation: 


<A,  “2  +i b 


where  Ax  represents  the  distance  between  points  on  this  array.  There¬ 
by,  Poisson's  equation  in  finite-difference  form  is  given  by 


(2-14) 


*1  "2  ^ 


=  Qc' 


(2-15) 
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where  Q  is  the  electrostatic  charge  and  is  given  by: 


o 


An  equation  similar  to  Eq.2-16  can  be  written  for  each  point 
on  our  approximating  grid .  When  this  is  done  we  have  1  00  poirts  and 

100  equations, with  100  unknowns - the  electrostatic  potential  at  each 

location  of  the  approximating  grid. 

Such  a  system  of  equations  can  be  written  in  matrix  form  yield¬ 
ing  a  tri-diagonal  matrix  of  rank  100.  There  is  a  multitude  of  different 
ways  one  can  solve  this  matrix  problem  ,  each  method  yielding  particu¬ 
lar  advantages  and  disadvantages. 

A  particular  source  of  difficulty  in  finite-difference  analysis  is 
round-off  error  when  solving  large  systems  of  equations .  For  this  rea¬ 
son  we  have  adopted  the  choleski  decomposition  technique[l4-16]|Which 
has  been  found  exceedingly  stable  in  the  presence  of  round-off  error. 

The  choleski  technique  is  based  upon  the  fact  that  if  a  matrix  B 
is  positive  definite,  than  B  can  be  written  as  the  product  of  a  lower  tri¬ 
angular  matrix  L  and  its  transpose  L*  this 

B  =  LL1  (2-17) 

and  we  have  a  two-step  solution  for  this  matrix  problem. 

It  is  emphasized  that  despite  a  direct  solution  for  this  matrix 
problem,  solving  Poisson's  equation  remains  an  iterative  process.  The 
exact  mobile  hole  distribution  is  unknown  and, therefore,  iteration  between 
the  solution  for  <Mx)  and  p(x)  must  continue  until  convergence  is  attained. 

In  an  attempt  to  accelerate  convergence  of  Poisson's  equation  we 
make  use  of  Gummels  algorithm  [^71  .  This  algorithm  is  based  upon  a  linear 
approximation  of  the  exponential  relations  appearing  in  the  Boltzmann  equa- 
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tion  for  holes,  Eq,2-2b.  We  assume  (/,m,  and  pm  represent  values  of 
electrostatic  potential  and  holes  respectively,  from  iteration  step. 

Further,  let  V  =  ' p  +  o  where  5  represents  a  small  quantity. 

Thus,  from  the  Boltzmann  relation 

(2-18) 

(2-19) 


exp(  V  m+i)  =  exp  {4>m)  £  1  +  8  m+*  J 


thus 

m+1  m  ,,  c  m+1 
p  =  n  (l-o 


representing  a  prediction  for  the  hole  density  at  any  location,  resulting 
from  the  next  iteration.  Thereby,  using  Poisson's  equation  we  have: 


a2  *  m+1 

dx^ 


K( 


.  m,  ,  m+1 
(p  )  <1> 


[n- 
“o  L  A 


m  m  .  m 
P  "P  * 


(2-20) 


In  finite  difference  form  only  the  center  terms  need  be  modified:  one 
term  added  to  the  finite  difference  element  and  a  second  term  to  the 
charge  element. 

In  calculations  of  electric  current  continuity  it  is  customary  to 
employ  similar  finite-difference  approximations  to  the  expression  for  elec¬ 
tric  current,  Eq.2-6,  and  substitute  the  results  into  an  expression  for 
the  divergence  of  electric  current.  However,  it  can  be  shown  that  this 
procedure  leads  to  numerical  instability  whenever  the  voltage  difference 
between  two  mesh  points  exceeds  about  5.1  x  10  ^  volts  [l8]  .  Following 
Sharfetter  and  Gummel,  the  expression  for  electric  current, 


I 


P 


-  q  p  P 

p 


dt// 

dx 


(2-21) 
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is  treated  as  a  differential  equation  in  p  with  J  ,  n  and  d  <A/dx  as¬ 
sumed  constant  between  mesh  points  . 

Thereby ,  we  have  a  differential  equation  of  the  form 


dp 

dx 


+  Ap  +  B  =  0 


(2-22) 


which  can  be  solved  between  two  finite-difference  mesh  points  yielding 


P  1 1  =  P„  exp  (U  -  U  ) 
n+1  n  n+1  n 


J  Ax 


"  (u„+i  -  u„» 


exp  (Un+1  -  U„M 


(2-23) 


where  U  represents  the  normalized  electrostatic  potential  (U  =  q  </»/  k  T). 
Equation  (2-23)  relates  the  hole  density  at  mesh  point  n+1  and  n  to  the 
electric  current  between  these  mesh  points.  From  Eq.(2-23)  we  can 
therefore  write 


1=  V 


p  n  AU 

n 


p  #x.\U 


Ax  j^exp(AU)  -  lj  Ax  ^exp(-AU)  -  1 


(2-24) 


where  AU  =  U  ,  -  U  .  Equation  (2-24)  expresses  the  electric  cur- 
n+1  n 

rent  between  mesh  points  n  and  n+1  in  terms  of  the  potential  and  hole 
densities  at  these  locations. 

Designating  the  electric  current  into  node  n  by  Jw  and  out  of  node 
n  by  J  ,  the  divergence  of  electric  current  at  node  n  is  given  by: 


dx 


w 


C 

n 
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where  C  represents  the  recombination/generation  term  assigned  to  node 
n.  After  substituting  Eq.2-24  into  Eq.2-25  for  the  east  and  west  elec¬ 
tric  currents  we  obtain  an  expression  for  the  divergence  of  electric  cur¬ 
rent  at  each  node  in  terms  of  the  hole  density  at  this  node,  and  the  two 
adjacent  nodes. 

Thereby  we  obtain  a  system  of  equations  for  the  hole  density - 

one  equation  for  each  node.  This  system  represents  a  matrix  equation 
for  electron  density  that  contains  a  sparce,  tri-diagonal  matrix.  This 
system  of  equations  is  solved  using  the  technique  of  Gauss  elimination; 
by  reducing  the  matrix  to  a  triangular  form  and,  thereafter,  solving  for 
the  electron  densities  through  backward  substitution. 

5 . 0  Boundary  Conditions 

Boundary  conditions  for  this  calculation  are  introduced  in  two  dif¬ 
ferent  arrays;  first,  in  the  electrostatic  potential  array  and,  second,  in 
the  hole  density  array.  Recognizing  that  all  essential  potential  drop  with¬ 
in  a  Schottky  barrier  device  is  encountered  within  the  semiconductor  mate¬ 
rial,  this  situation  is  introduced  through  an  adjustment  of  the  boundary  con¬ 
ditions.  Additionally,  this  mathematical  model  is  directed  toward  an  ef¬ 
fective  generation  of  holes  at  the  metal-semiconductor  interface;  such  a 
mechanism  is  also  introduced  through  an  adjustment  of  boundary  condi¬ 
tions  at  this  location. 

For  the  solution  of  Poisson's  equation  we  assume  both  end  points 
of  the  one-dimensional  array  are  maintained  at  a  specified  electrostatic 
potential:  zero  at  the  barrier  interface,  and  at  some  positive  applied 
voltage  at  the  semiconductor  ohmic  contact.  At  each  of  these  locations 
we  have  a  finite-difference  equation  of  the  type: 

A  +  B  +  D  t//  =0  (2-26) 

w  c  e  c 
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If  we  designate  the  barrier  interface  at  the  western  most  location  and 
the  ohmic  contact  as  the  eastern  most  location,  it  is  clear  that  A  at 
the  barrier  and  D  at  the  ohmic  contact  do  not,  in  fact,  exist.  Thus, 
setting 

(Barrier  boundary)  D  =  O,  B-=  1,  Q  =  (-N  +  p) 

C  A 

(ohmic  contact)  A  =  O.B  =  l,Q  =  V 

c  appl 

The  boundary  conditions  for  the  potential  array  are  readily  satisfied. 

In  the  presence  of  an  electric  current  the  hole  quasi-Fermi  po¬ 
tential  has  some  unspecified  gradient  distribution  between  the  ohmic 
contact  and  the  barrier  interface.  This  gradient  distribution  is  unknown, 
and  it  will  depend  upon  the  electric  current  density  at  each  location  with¬ 
in  the  device.  If,  indeed,  the  semiconductor  ohmic  contact  lies  far  re¬ 
moved  from  the  diode  space-charge  layer,  the  quasi-Fermi  potential  at 
that  location  will  equal  Fermi-potential,  and  at  a  point  that  is  charge 
neutral.  Further,  if  the  hole  Fermi-potential  is  maintained  constant  from 
the  ohmic  contact  to  the  Schottky  barrier  we  would  model  a  device  of  this 
type  .n  the  absence  of  an  electric  current.  Thus,  reasonable  boundary 
conditions  appear  to  be  attained  by  maintaining  at  the  ohmic  contact  a 
hole  quasi-Fermi  potential  equal  to  the  Fermi- potential  for  the  p-type 
semiconductor  material.  At  the  barrier  interface  we  calculate  the  hole 
density  using  Eq.1-18,  and  thereby  establish  the  quasi-Fermi  potential 
at  that  location. 

Thus,  with  knowledge  of  the  quasi-Fermi  potentials  at  the  ohmic 
contact  and  at  the  barrier  interface  we  can  readily  calculate  the  electron 
densities  at  these  locations  and,  thereby,  establish  boundary  conditions 
for  solving  the  equation  of  electric  current  continuity . 
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6 . 0  Iterative  Solution  for  the  Schottky  Barrier  Diode 

As  previously  discussed,  iterative  techniques  are  used  to  solve 
this  problem  of  the  Schottky  barrier  diode.  A  given  solution  is  initiated 
by  a  first  "guess"  concerning  the  electrostatic  potential  and  the  hole 
density  distributions .  Thereafter,  using,  iterative  techniques ,  Poisson's 
equation  is  solved ,  and  after  each  solution  the  hole  distribution  is  up¬ 
dated,  maintaining  the  initial  guess  for  the  hole  quasi-Fermi  potential; 
by  recalculating  the  hole  density  we  are  establishing  an  improved  solu¬ 
tion  for  the  electrostatic  charge  distribution. 

Prior  to  each  new  calculation  of  the  electrostatic  charge  in 
Poisson's  equation,  a  determination  is  made  concerning  the  degree  of 
convergence.  This  is  accomplished  by  calculating  for  each  location  of 
the  potential  array  the  expression: 


A  0  +  B  0  +  D0  -  Q  =  R 
e  c  w  c 


(2-27) 


which  yields  the  error,  or  residual,  at  a  given  spatial  location.  Ideally, 
the  true  solution  for  this  problem  would  yield  a  residual  of  zero  everywhere. 

Termination  of  this  particular  aspect  of  the  calculation  is  attained 

-4 

when  the  maximum  residual  upon  entry  to  the  potential  calculation  is  10 
(or  less).  Experience  shows  that  any  attempt  at  further  refining  the  solu¬ 
tion  is  of  little  practical  value.  In  practice,  this  initial  convergence  of 
Poisson's  equation  is  obtained  with  100  (or  less)  iterations  between  elec¬ 
trostatic  potential  calculation  and  the  hole  density,  for  a  given  quasi- 
Fermi  potential  calculation. 

After  attaining  convergence  for  this  initial  part  of  the  potential 
calculation  we  obtain  a  very  accurate  solution  for  Poisson's  equation, 
for  our  first  guess  about  the  quasi-Fermi  potential  distribution.  Next, 
an  improved  solution  for  the  quasi-Fermi  potential  is  attained  by  solving 
the  continuity  equation  for  the  hole  density  and  recalculating  a  new  quasi- 
Fermi  potential. 
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After  the  initial  convergence  of  Poisson's  equation  we  no  longer 
require  the  iteration  between  the  electrostatic  potential  calculation  and 
the  hole  density,  maintaining  the  quasi-Fermi  potential  constant.  After 


each  new  solution  of  Poisson's  equation  we  immediately  recalculate  a 
new  hole  distribution  and,  hence,  a  new  hole  quasi-Fermi  potential. 

Final  convergence  for  this  sytem  of  equations  is  based  upon 
simultaneously  attaining  thFee  differnt  convergence  criterion: 


1.  When  the  starting  residual  is  low  prior  to 
entry  into  the  electrostatic  potential  cal¬ 
culation. 

2 .  When  the  starting  residual  is  low  prior  to 
entry  into  the  calculation  of  electric  cur¬ 
rent  continuity. 

3 .  When  electric  current  continuity  is  attained 
within  1%  of  the  total  electric  current. 


From  a  practical  point  of  view,  excellent  convergence  is  usually  at¬ 
tained  within  100  iterations  between  the  solution  of  Poisson's  equation 
and  the  solution  for  electric  current  continuity. 
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APPENDIX  1 


Computer  Programs  Developed  E)uring  This  Investigation 


Program  Function  Page 

A.  Calculating  Fermi  Level  vs.  Temperature .  83 

B.  Calculating  Unoccupied  States  Between  E  and  E. .  9 7 

c  I 

C.  Calculating  Unoccupied  States  Between  E.  and  E. .  115 

b  f 

D.  Calculating  Unoccupied  States  Between  E  and  E, .  133 

c  b 

E.  Calculating  Transition  Energy  Ec .  151 

F.  Calculating  Saturation  Currents,  Jt,  J  ,  and  J .  171 

12c 

G.  Calculating  the  Volt-Ampere  Characteristics  of  a 

Schottky  Barrier .  195 

H.  Calculating  Many  Equilibrium  Properties  of  Silicon .  223 

I.  Numerically  Solving  Poisson's  Equation  . .  229 


Note:  All  computer  programs  listed  in  APPENDIX  1  are  self  contained. 
The  listings  within  a  given  function  group  can  be  keypunched , 
compiled,  linked,  and  run  as  a  unit. 
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APPENDIX  1-A 


Program  for  Calculating  the  Fermi-Level  vs.  Temperature 


Title  Page 

MAIN .  84 

EFERMI .  86 

ROOT .  88 
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XINT 
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IF(XX.6T.1.fD2)  GOTO  30 

CALL  OUAD1 (X INI, XLBND.XUBND, RES, 14, DUMMY) 
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Fig. 31  Ionized  Impurity  Atom  (Boron)  Density  in  Silicon 
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